The Higgs mode in a superfluid of Dirac fermions 
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We study the Higgs amplitude mode in the s-wave superfluid state on the honeycomb lattice 
inspired by recent cold atom experiments. We consider the attractive Hubbard model and focus on 
the vicinity of a quantum phase transition between semi-metal and superfluid phases. On either side 
of the transition, we find collective mode excitations that are stable against decay into quasiparticle- 
pairs. In the semi-metal phase, the collective modes have "Cooperon" and exciton character. These 
modes smoothly evolve across the quantum phase transition, and become the Anderson-Bogoliubov 
mode and the Higgs mode of the superfluid phase. The collective modes are accommodated within 
a window in the quasiparticle-pair continuum, which arises as a consequence of the linear Dirac 
dispersion on the honeycomb lattice, and allows for sharp collective excitations. Bragg scattering 
can be used to measure these excitations in cold atom experiments, providing a rare example wherein 
collective modes can be tracked across a quantum phase transition. 
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Introduction- Spontaneous symmetry breaking of con- 
tinuous symmetries gives rise to two typical collective 
excitations - gaplcss Goldstone modes and a gapped am- 
plitude mode, also called the Higgs mode [l|. While the 
Goldstone mode has been observed in various contexts, 
the Higgs mode has evaded observation with rare excep- 
tions such as NbSe2, which has coexisting charge den- 
sity wave and superconducting order 0, Q and multifer- 
roic Ba2CoGe207 Remarkably, two recent experi- 
ments have successfully observed this mode by tracking 
collective excitations across a quantum phase transition. 
The first involves pressure studies of TICUCI3, a mag- 
netic material which undergoes a transition from dimcr 
order to magnetic order [f|. The second is the realiza- 
tion of the Bose-Hubbard model in ultracold gases, with 
a visible amplitude mode near the superfluid-Mott tran- 
sition [H, 0- In this letter, we propose a novel scheme 
to observe the Higgs mode in a Fermi superfluid. Hith- 
erto, the Higgs mode has never been observed in Fermi 
supcrfluids as it decays into pairs of quasiparticles. Our 
proposal circumvents this issue by exploiting a special 
feature of the honeycomb lattice geometry which allows 
for a window in the quasiparticle-pair continuum - the 
Higgs mode survives as a stable excitation inside this 
window. 

Inspired by the recent realization of the honeycomb 
optical lattices in cold atom experiments @, we study 
the attractive Hubbard model in this geometry: 
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Parameter tij denotes hopping amplitude between near- 
est neighbor (tij — t) and next-nearest neighbor sites 
(tij = t'). U is an on-site attractive interaction and 
li is the chemical potential. We envisage a setup with 
a deep optical lattice to trap two hyperfine species of 
fermions, and a magnetic field on the attractive side of 



a Feshbach resonance @. This model hosts a superfluid 
state of Dirac fermions, with several interesting implica- 
tions [HI E| ■ In this proposal, we make use of two key 
features: (i) strictly at half-filling, there is an interaction- 
tuned quantum phase transition from a semi-metal phase 
to an s-wave superfluid. This has been demonstrated by 
sophisticated quantum Monte Carlo simulations on very 
large system sizes [lU, EH- This transition is a conse- 
quence of the Dirac cone dispersion which leads to vanish- 
ing density of states at the Fermi level, thereby necessi- 
tating a critical interaction strength to induce superfluid 
order [Uj]. (ii) In the semi-metal phase, the two- 
particle continuum has a window structure, again a con- 
sequence of the Dirac cone dispersion [l5| ■ A collective 
mode excitation propagating inside this window is stable 
against decay into quasiparticle-pairs. We show that this 
window structure persists in the superfluid phase, thus 
allowing for a stable Higgs mode excitation. 

Fig. [TJb) shows the phase diagram of this model at 
half- filling. Our key findings are summarized in Fig.[TJc)- 
(f): (i) there are sharp collective mode excitations on ei- 
ther side of the transition. The two-particle continuum is 
shown as the shaded region, note the window structure. 

(ii) In the semi-metal phase, there are three degenerate 
collective modes with "Cooperon" and exciton character. 

(iii) On the superfluid side, there is a Goldstone mode 
and remarkably, a distinct superfluid amplitude (Higgs) 
mode. These modes can be observed in a cold atoms 
experiment using Bragg scattering. This is a rare exam- 
ple wherein relevant collective excitations can be tracked 
across a quantum phase transition. 

Mean field theory- The ETH group Q has studied 
fermions loaded onto a honeycomb optical lattice with 
tunable anisotropy. We consider the attractive Hub- 
bard model in the isotropic honeycomb lattice. As dis- 
cussed in Ref. [ll| , the isotropic limit is expected to have 
the highest superfluid transition temperature and is the 
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FIG. 1: Isotropic honeycomb lattice with basis vectors (a). 
Phase diagram of the attractive Hubbard model on a honey- 
comb lattice at half-filling with t' = as obtained from the 
mean-field theory (b). Evolution of the elementary excita- 
tions across the quantum critical point (QCP) of the semi- 
metal (SM) to superfluid (SF) phase transition (c)-(f). The 
dash-dot line in each panel and the dash-dot-dot lines in (f) 
show the asymptotic dispersions of the continuum edge and 
the superfluid collective modes for small q, respectively. 



most promising for experimental realization. We decom- 
pose the Hubbard interaction in the superfluid channel 
using the order parameter A = U(ciiCif), taken to be 
real. For brevity, we introduce a vector operator con- 
sisting of creation and annihilation operators &(p) = 

(cp,a,t) c Lp,a,4.' c p,f>,t' c -p.6,i)* ( a an d b denote the two 
sublattices as shown in Fig. [Ha)). The mean field Hamil- 
tonian can be written as Hmf = J2 P ^Hp)h(p)^(p)> 



where h(p) = {x p + Re(7p)<7 x - Im(7p)cr y }T 3 - An, 
7p = -i(l + e lp ai + e lp - a2 ) and x p = -2£'[cos(p • ai) + 
cos(p • a 2 ) + cos(p ■ (ai — a 2 ))} — p, with a\ and a 2 being 
the two basis vectors shown in Fig. Q] (a). We take the 
lattice spacing to be unity, f and a are the Pauli ma- 
trices in the Nambu and sublatticc space, respectively. 
The single-particle Green's function for the mean-field 
Hamiltonian is given by 
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Here, we denote p = {p, iw„), where w„ is the fermion 
Matsubara frequency. The gap and number equations 
are obtained from the off-diagonal and diagonal elements 
of the Green's function G vv as (hereafter, we restrict 



ourselves to zero temperature) 
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where E ± (p) = y (£p ) 2 + A 2 is the spectrum of the Bo- 

goliubov quasiparticles, £p = x p ± |7 P |, and TV is the 
number of lattice sites. At half-filling, the self-consistent 
solution of A becomes non-zero for U > U c indicating a 
transition from semi- metal to superfluid phases [ToL fill ] . 
For t' = 0, mean-field theory gives U c ~ 2.23i. Quantum 
Monte Carlo gives the same transition, except with U c 
renormalized to ~ 3.869 [l2| . In the rest of this letter, 
we use mean-field results with the understanding that 
fluctuations will renormalizc U quantitatively. We note 
that U c weakly depends on the value of t' . 

Generalized Random Phase Approximation ( GRPA )— 
On either side of the critical point, there are low-lying 
density and pairing fluctuations. We use a generalized 
random phase approximation (GRPA) scheme to evalu- 
ate density and pairing response functions. We follow the 
Green's function approach of Cote and Griffin We 
denote matrix susceptibilities containing the response to 
weak density and pairing perturbations, respectively, as 
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where q = (q, iO„) (f2„ is a boson Matsubara frequency). 
Any susceptibility \ is defined as 



-n„ri 2 ) 



X/T(«) = -E / dr 12 (T T 5f(l)8g(2))e~^ 

(7) 

where 1 = (ti 1 ,Vi,t{) (h denotes the unit cell, v\ the 
sublattice, and t\ an imaginary time), T\ 2 = ri t — r; 2 , 
T12 = Ti — t 2 and Sf = f — (/). The density and pair 
annihilation operators are written as n = n-\ + n\. and 
m = c±Cf, respectively. The GRPA equations read [Tl], 



IT J 

A^(q) = A°^(q) + — J2J2^ U3 (j> + q) 

v 3 p,U n 

xA U3U2 (q)G U3Vl (p), (8) 
A v ^(q) =A v ^{q)~U^L^ i {q) r lT{A^-{q)}, (9) 



where A is either L or M and G^Hp) = T 3 G V1U2 (p) . 
A denotes the bare susceptibility [17j . A includes an 
infinite sum over ladder diagrams, while A is the final 
result which also includes bubble diagrams. 
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Undamped Higgs mode- In the superfluid phase, we 
solve GRPA equations (fS) and © to evaluate the 
amplitude and phase correlation functions Xaa 2 (?) = 

1V2 M+xZZM) and Xe T 2 (?) = &(x2Z(q)~ 

The amplitude and phase fluctuation op- 
erators are given by 6 A = ^{8m + 5mt) and 86 = 
2^(i5m — Sm^), respectively. For the case of t' = 0, 
the expressions simplify and we can identify their respec- 
tive poles, which we denote "Higgs" and "AB/Leggett" . 
These poles satisfy 



Higgs: 1 = -{G + D) + \R\, 



(10) 



AB/Leggett: — = -{G - D) + ^AF 2 + \R\ 2 . (11) 



We have defined 



c = -Y 



E + E' 



N ^ (iVt n ) 2 - {E + E') 2 ' 



D = -Iv) * 
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F = — V A 



N ^ E {ifl n ) 2 -(E + E 
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(12) 
(13) 



N ^ E'E (in n ) 2 — (E + E') 2 ' 



N ^ E'E {in n ) 2 -{E + E') 2 ' 



(15) 



Here, we have denoted E = E(p), E' = E(p+q), 7 = j p , 
and 7' = 7 p +q. On the other hand, solving Eqs. ^ and 
([9]) for density response, we find that Xee" 2 (q) ^ Xnn 2 (<?) 
when i' = 0. Thus, the density response function only 
retains the AB/Leggett pole given in Eq. (fTTj) . 

Setting q = Q n = in Eq. (fTT|) . we recover the gap 
equation ([3]). Thus, the superfluid phase has gapless col- 
lective mode(s) arising from phase fluctuations. In fact, 
at half- filling, the AB/Leggett pole in Eq. (fTTj) is a double 
pole corresponding to two gapless modes: the Anderson- 
Bogoliubov (AB) mode and the Lcggett mode [ll|. The 
AB mode is the usual Goldstone mode associated with 
U(l) symmetry breaking [18L Il9j. The Leggett mode is 
composed of out-of-phase fluctuations between sublat- 
tices [ITj] - it acquires a gap away from half- filling fllj ) . 
The AB and Leggett modes become degenerate at half- 
filling reflecting a special pseudospin SU(2) symmetry of 
the Hubbard model [11]]. For small q <C 1, Eq. (fTTj) gives 



the dispersion relation of the AB/Leggett mode to be [2lJ 



ujab = Au_Fg, A 
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(16) 



be understood using the mechanical analog of motion 
along the radial direction of the famous "Mexican hat" 
potential; the energy gap stems from the finite curvature 
of the potential along the radial direction. 

Remarkably, the Higgs mode disperses below the quasi- 
particle pair continuum in Figs, flje) and (f). In partic- 
ular, close to the M point, it is well separated from the 
lower edge of the continuum. This is to be contrasted 
with the case of typical supcrfluids: due to the under- 
lying Fermi surface, the continuum exhibits a horizontal 
edge near q ~ f3(. The Higgs mode therefore enters 
the continuum, becomes heavily damped and is unob- 
servable. 

In our case, the Higgs mode in Fig. [JJ is undamped 
over large sections of the Brillouin zone. For q <C 1, solv- 
ing Eq. (|10[) . the Higgs mode has the dispersion relation 
w Hi gg s = 4A 2 + v 2 F q 2 [2l|. The Higgs mode disperses 
below the lower edge of the continuum which is given by 
Wedge = mm p [(E(p)+E(p+q)) 2 } ~ 4A 2 + 2v 2 F q 2 (q « 1). 
This window or arch in the continuum, shown in Fig. [TJ is 
a consequence of the Dirac-likc dispersion of underlying 
fcrmions. The Higgs mode stays undamped as long as it 
lies within this window. Even if we go slightly away from 
half-filling, the Higgs mode survives undamped. Close to 
?i ~ 0.9 oin~ 1.1, the window disappears because of the 
presence of the Fermi surface - as a result, the Higgs mode 
is strongly damped. The Higgs mode, the AB mode, and 
the lower edge of the continuum become degenerate at 
the QCP for q < 1: uj A b = ^Hi gg s = w odgo = v F q. 

The AB mode and Leggett mode are strongly coupled 
with density fluctuations; they appear as poles in the den- 
sity response function Xnn" 2 GO- However, when t' = 0, 
the Higgs mode has no corresponding pole in the density 
response function. Thus, the Higgs mode is composed of 
pure amplitude fluctuations and cannot be excited by a 
density perturbation. This reflects the underlying SU(2) 
pseudospin symmetry [ll| in the problem. A small finite 
t' breaks this symmetry and forces the Higgs mode to 
acquire a density component: the density response then 
shows a peak at the Higgs mode, as shown in Figf2j 

Collective modes in the semi-metal phase: Cooper- 
ons and excitons- In the semi-metal phase, set- 
ting A = 0, density and pair response func- 
tions become decoupled in the GRPA equations (f5j)- 
©. The susceptibilities satisfy the usual RPA equa- 
tions x VlV2 M = x^Hq) + UJ2„ x 0viV ?{q)x V3V2 M, 
xZ^(q) = xZT(q) - uY.„ 3 xZT(q)xT 2 Aq)- The 

bare susceptibility X° mm i describes a single rung diagram 
with particle-particle (hole-hole) excitations, and Xn n 
describes a single bubble diagram with particle-hole ex- 
citations. They are given by 



where vf = 3i/2 is the Fermi velocity at the Dirac points. 

Setting (q,iVL n ) = (0, 2A), Eq. (fTUj) also reduces to the 
gap equation. Thus, there exists a gapped collective mode 
with the energy gap 2A at q = 0. This is the 'Higgs' 
mode or the amplitude mode Q arising from amplitude 
fluctuations of the superfluid order parameter. It can 



X°:,L7(q) 



An^n 



(q) 



where ftp 1 " 2 



1 






^"X-p 


2N 






1 








27V 


£ + , -£~-i£7„ 


-fj— iOn 

^>p + q ^P (L 



4- P^vS u 4- p-tyvfi ufi 



,(17) 
,(18) 
and 



4 



"J/1 O u i 



We denote 



^ = 7 P /l7pl 



From the denominator in Eq. (|17[) , we see that pairing 
response arises from particle-particle (or hole-hole) ex- 
citations. An undamped pairing mode, occurring below 
the particle-particle continuum in Fig. Q3c), is therefore 
a two particle bound state with well defined momentum 
and energy. It can be understood as a preformed Cooper 
pair - we call this a Cooperon excitation (22j . Similarly, 
from Eq. (fT5|) . we see that density response arises from 
particle-hole excitations. An undamped density mode is 
thus a particle-hole bound state - we call this an exciton. 

The dispersions of Cooperons and excitons are deter- 
mined by the poles of the corresponding response func- 
tions, giving \i-Ux mm i(q)\ = and |J + Cx°„„(g)l = °- 
With t' = 0, these reduce to the identical equation 



(1 + Ua(q)) 2 - 



U 2 \f3(q)\ 2 
|7pI + 



= 0, 

\lp+q\ 



n an 



2 _ 



(|7 P I + l7p+«l) 
■*- ) (|7pI + |7p+,I) 
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l7p+ql 



(19) 
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Thus, the Cooperon and the exciton are degenerate when 
t' = 0. Their dispersion is shown in Figs. Q] (c) and (d) - 
the modes are undamped as they lie below the two par- 
ticle continuum. In particular, they are well separated 
from the continuum in the vicinity of the M points. We 
suggest that experiments should probe this region to ob- 
serve the collective excitations. This feature of the M 
points can be understood from the single particle band 
structure which has saddle points at these wavevectors. 
They consequently have a very large density of states 
which provides large phase space for the Hubbard inter- 
action to form two-particle bound states. 

These collective modes in the semi-metal phase were 
predicted many years ago - using an insightful single- 
cone approximation - Ref. 15| reported a triplet exciton 



mode in the repulsive Hubbard model. The authors iden- 
tified the window structure in the continuum as capable 
of accommodating stable modes. Mapping their results 
to the attractive Hubbard case [l3|, the triplet excitons 
translate to Cooperon and exciton modes. We reaffirm 
their prediction starting from a microscopic picture tak- 
ing into account the sublattice structure. Our expressions 
also agree with those of Ref. [23| - which only considers 
the r — K segment and concludes that there is no un- 
damped mode. However, we find an undamped mode in 
the r — M and M — K directions. 

Cooperon condensation- As we approach the crit- 
ical point from the semi-metal side, the energy of 
the Cooperon and exciton decreases progressively (see 
Fig- IUc) and (d)). Precisely at the transition, the 
Cooperon "softens" at q = and undergoes condensa- 
tion. In fact, setting q = fi n = 0, the Cooperon pole in 
Eq. (|T9"|) reduces to the gap equation Since Cooper- 
ons and excitons are degenerate for t' = 0, the exciton 
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FIG. 2: Intensity of dynamic structure factor correspond- 
ing to density response S(q,u>) = —lm[xnn(q,u})]/'K for t' = 
—0.05t (a). The cross section for the momentum at the M 
point (b). The upper peak corresponds to the Higgs mode. 



can also condense at the critical point. That gives rise to 
the sublattice-CDW state - which is degenerate with the 
superfluid state due to SU(2) pseudospin symmetry. For 
t' 7^ 0, this degeneracy is lifted in favour of the superfluid 
and the Cooperon condenses preferentially. 

As we cross U c and enter the superfluid phase, Cooper- 
ons and excitons hybridize to become the AB, Leggett 
and Higgs modes. The excitonic component, when 
present, allows these modes to have peaks in the density 
response function. The Cooperonic component manifests 
as peaks in the pairing response. Thus, the collective 
modes evolve smoothly across the QCP and carry signa- 
tures of the underlying spontaneous symmetry breaking. 

Visibility of the Higgs mode- Our Higgs mode is stable 
against decay into pairs of fermions due to the window 
structure in the two particle continuum. However, if we 
go beyond RPA level, the Higgs mode may decay by emit- 
ting AB/Leggett excitations which are lower in energy. 
We argue that this merely leads to broadening of the 
Higgs mode. In our semi-metal to superfluid transition, 
due to the pseudospin symmetry present when t' = 0, 
the order parameter can be thought of as an 0(3) object. 
Quantum Monte Carlo simulations of the 0(3) ordering 
transition show that the Higgs mode survives although 
it is broadened by decay into Goldstone bosons [24[ . We 
expect this to be true in our case as well. 

We suggest Bragg spectroscopy measurements on a 
Fermi superfluid in a honeycomb optical lattice as a 
way to measure this undamped Higgs mode. In this 
technique, a two-photon process imparts a density- 
"kick" to the system. The response to this pertur- 
bation can be quantified by measuring the momentum 
transferred or the energy absorbed. The momentum 
transferred is a measure of the dynamic structure fac- 
tor related to the density response function S(q, uS) = 
— Im[x„„(<7, uj)]/tt [2|| - it can detect collective modes as 
long as they have a density component. For any small 
t' ^ 0, the Higgs mode in our model has a density com- 
ponent which makes it visible to Bragg spectroscopy. A 
small t hopping is expected to be present in an optical 
lattice setup any way (2(|. Figure [2] shows S(q,w) for 
t' = —0.05. The sharp intensity peak for the Higgs mode 
can be clearly seen below the continuum. An alternative 
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approach is to measure the energy absorption in response 
to a weak shaking of the optical lattice [7| . 
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I. SUPPLEMENTARY MATERIALS 
A. Bare susceptibility in GRPA 

To solve the GRPA equations (8) and (9), the bare 
susceptibilities L° and M° are evaluated using the mean- 
field Green's function in Eq. (2) to give 
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where we introduced the tensor 
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Following Rcf. |16| . we introduce a column vector A{- 
C,A4) and a 4 x 4 matrix I? as 
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The GRPA equations are cast into the form 

A™ (q) = A " 1 " 2 (q) + Uj2 'D VlV3 {q)A V3V2 (q)(27) 
A VlU2 {q)=A VlV2 (q) - U^2C VlV3 (q)A V3V2 {q). (28) 



(24) The above equations are easily solved to give 



A aa {q) A ab (q) 
A ba (q) A bb (q) 

A ab (q) 
A bb (q) 



A. aa (q) 
A ba (q) 



1 + UL aa {q) 
UL ba (q) 

i-UV aa (q) 
-UV ba {q) 1 



UL ab _(q) 
1 + UL bb {q) 

-UV ab (q) 
- UV bb (q) 



A_ aa (q) A_ ab {q) 
A ba (q) A bb {q) 

( A 0aa (q) A 0ab (q) 
{A oba (q) A obb (q) 



(29) 
(30) 



r 



Here, we denoted A h 



Tr{A Ul1 ' 2 } = A\ 



2 +A V 4 1V2 . 



Here, we have replaced iCl n — >• uj. At q = 0, Eq. ([32 
reduces to 



B. Energy dispersion of the Higgs mode 

We derive the analytic expression of the energy dis- 
persion of the Higgs mode for small momentum ij < 1 
following the approach of Ref . Q • The gap equation ([3]) 
can be rewritten as 



U 
2N 



i 

e 1 



= o. 



(31) 



We denote E(p) = E, E(p + q) = E' , 7 = j p , and 7' = 
7p +q . Subtracting Eq. flST]) from 1 + U[(C+D)-\R\] = 
which is equivalent to Eq. ([TO]) , we obtain 



U(X-\R\) = 0, 



(32) 



where 



X= (C + D) 



N ^ 



2N 



E + E 1 
EE' 



— V 

2N ^ 



E + E 1 uj 2 - (| 7 | 2 + | 7 '| 2 ) -4A 2 



EE' 



uj 2 -{E + E') 



l\2 



(33) 



v 1 1 



N ^ Euj 2 -4E 2 
v 

r-Zt 

= V I de 



P{e) 







y/e* + A 2 [w 2 - 4(e 2 + A 2 ) 



0, (34) 



where y = uj 2 /A- A 2 and p{e) = jjJ2 P S { £ ~ W) is the 
density of states of the fermion energy band. If we set 
uj = 2A, the denominator in the integrand of Eq. (|34|) is 
proportional to e 2 , while the numerator is proportional 
to e for small e because p(e) oc e. The integral in Eq. (f34|) 
is thus well defined in the limit uj — > 2 A. In this limit, 
Eq. (f34"l) is satisfied when y = and consequently the 
energy of the Higgs mode is obtained as WHiggs(q = 0) = 
2A. 

To derive the energy dispersion for small q, we expand 
Eq. ([32]) to second order in q. Using the relations 
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lp+q ~ 7p + <5 7l + J 72 , <5 7l = -it{e ip ' ai (Q ■ 01) + e ip ° 2 (<Z ■ a 2 )}, (35) 

^72 = ^{e lp ai (g • ai) 2 + e Jp ° 2 (g • a 2 ) 2 }, (36) 

wi = Re[7*c>7i] = -t 2 {(sin(p • a\) + sin(p • a 3 ))(g • a\) + (sin(p • a 2 ) - sin(p • a 3 ))(g • a 2 )} , (37) 

= Im[7*57i] = t 2 {(1 + cos(p • ai) + cos(p • a 3 ))(q ■ a x ) + (1 + cos(p • a 2 ) + cos(p • as))(g • a 2 )} , (38) 



u) 2 = Rc[7*#7 2 ] 



{(1 + cos(p • ai) + cos(p • a 3 ))(q • a x ) 2 + (1 + cos(p ■ a 2 ) + cos(p • a 3 ))(g • a 2 ) 2 } , (39) 



w' 2 = Im[7*(57 2 ] = -— {(sin(p • a x ) + sin(p • a 3 ))(g • ai) 2 + (sin(p • a 2 ) - sin(p • a 3 ))(g • a 2 ) 2 } 

l7 P +g| 2 — |7 P | 2 + si + s 2, si = 2u>i, s 2 = |<57i| 2 + 2w 2 , 
l^7i| 2 = t 2 {{q ■ ai) 2 + 2cos(p • a 3 )(g ■ ai)(g • a 2 ) + (g • a 2 ) 2 } , 

I 



(40) 

(41) 
(42) 



one finds that Eq. (|33|) becomes 

2 w 2 - (2| 7 | 2 + Sl + s 2 )-4A 2 



— V 

2N ^ E 
p 



1 w 2 



TV ^ E 
p 



UJ 2 - 4£ 2 

(2| 7 | 2 + S2 )-4A 2 



AE 2 



(43) 



Since the factor sin(p • a^) is odd for p. the summation 
for p including this factor vanishes. Similarly, R can be 
approximated as 



ReR 



j_ 2 ( 7 + <J 7l + <5 72 ) 7 * 



N ^ E 
p 

N ^ 



UJ 2 - 4£ 2 

2 |7| 2 + ^ 2 



(44) 



N ^ Elo 2 -AE 2 
p 



ImR = - 



1 E 



AT ^ E lo 2 - 4£ 2 
p 



(46) 



In evaluating further Eqs. (|4"3")l , (|4"5"j) . and (|4^|) . we en- 
counter the factor 



E 



1 COs(p • Oj) 

£ iJ 2 - AE 2 



(47) 



Here, we replace cos(p-eii) in the integrand with its value 
at the K (K 1 ) point, i.e., (cos(p ■ a{)) = cos(pk • ai) = 
— 1/2. At half-filling, since the Fermi level is at the K 
{K') point, this replacement is justified for small q. As a 
result, Ri and R 2 vanish and we finally obtain 



X-\R\^(Ay-vW)±Y,l 



N ^ Eu 2 -AE 2 
p 



0, (48) 



where vp = it/ 2 is the Fermi velocity. Consequently, the 
dispersion of the Higgs mode is obtained as 



<4i ggs -4A 2 + z; 2 g 2 . 



(49) 



In the limit A 3> vpq , WHi ggs is approximated 



2A + v 2 F q 2 /AA which is plotted in Fig. 1 (f). 



On the other hand, at the transition point with A = 0, 
the dispersion for small q coincides with that of the lower 
edge of the continuum as i^Higgs = vpq. 



C. Energy dispersion of the AB/Leggett mode 

We derive the energy dispersion of the AB/Lcggctt 
mode for small momentum. Subtracting Eq. pip from 
1 + U[(C — D) — ^JAF 2 + |i?| 2 ] = which is equivalent 
to Eq. (11), we obtain 



U 



(Y - ^JAF 2 + |i?| 2 ) = 0, 



Ro + R 2} (45) where 



E + E' 



2N EE' 
p 

E + E' ^ 2 -(| 7 | 2 + | 7 f ) 
V ^ EE' lj 2 -{E + E') 2 ' 



-E 

ON L^i 



Setting q = 0, we obtain 

1 uj 2 -2| 7 | 2 



Y 



1 E 



N ^ E uj 2 - AE 2 
p 



R = — 



F 



2 ^ 1 

i 

L E 



17 1 



AT ^ B w 2 - AE 2 
p 

1 Aw 



diui — 2d 2 , 



= -2d 2 



N ^ Elo 2 - AE 2 
p 



d 2 Au, 



where 



i i 



d2 M = ^E 



AT ^ E uj 2 - AE 2 ' 
p 

1 l7| 2 

N ^ E uj 2 - E 2 ' 
p 



(50) 



(51) 



(52) 
(53) 
(54) 

(55) 
(56) 



Thus, Eq. (JSDl) reduces to 

Y- V4F 2 + |i?| 2 

d 1 w a (d 1 w 2 -4d a -4d 1 A a ) 

= , = 0. (57) 

(diw 2 - 2d 2 ) + A 2 w 2 + 4d 2 

From d i= i 2 (to < 2 A) < 0, we obtain the gapless 
AB/Leggett mode: u>ab(q = 0) = 0. Note that the term 
within parentheses in the numerator of Eq. (|57[) at u) = 
is found to give 

- M 2 (lo = 0) - 4d 1 (w = 0)A 2 = I £ I = I. (58) 

p 

For small q, expanding i* 1 , i?, and Y to second order 
in q, we obtain 

Y ~ diw 2 - 2d 2 - diu|g 2 , (59) 
F-Adxo;, i?~-2d 2 . (60) 

In Eqs. (|59p and (|60[) . we used the same approximation 
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as the one for Eq. (|48|) . As a result, we obtain 



Y- v/4F 2 + |i?| 2 

= 0)(w 2 /[/ + 4d 2 (w = 0)u|(7 2 ) 
~ (diw 2 - 2d 2 - 2d 1 v 2 F q 1 ) + v /4A 2 d 2 w 2 + 4d 2 
= 0. (61) 

We have used Eq. (|58|) to derive the final expression. The 
pole is thus given by 

wab = A«f<7, (62) 

A 2 =4^(^ = 0)1 = ^^. (63) 

p 

Note that A < 1 from the gap equation (2). At the 
transition point (A = 0), A = 1 and thus the AB/Leggett 
mode becomes degenerate with the Higgs mode as well 
as the edge of the continuum as wab = vfQ- 



